Exchange bias phenomenology and models of core/shell nanoparticles 
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Some of the main experimental observations related to the occurrence of exchange bias in mag- 
netic systems are reviewed, focusing the attention on the peculiar phenomenology associated to 
nanoparticles with core/shell structure as compared to thin film bilayers. The main open questions 
posed by the experimental observations are presented and contrasted to existing theories and models 
for exchange bias formulated up to date. We also present results of simulations based on a simple 
model of a core/shell nanoparticle in which the values of microscopic parameters such as anisotropy 
and exchange constants can be tuned in the core, shell and at the interfacial regions, offering new 
insight on the microscopic origin of the experimental phenomenology. A detailed study of the of the 
magnetic order of the interfacial spins shows compelling evidence that most of the experimentally 
observed effects can be qualitatively accounted within the context of this model and allows also to 
quantify the magnitude of the loop shifts with striking agreement with the macroscopic observed 
values. 
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I. INTRODUCTION 

Magnetic fine particles have attracted a constant 
interest among the scientific community during the 
last decades because of their increasing number of 
applications^. The demand for miniaturization and 
the availability of new synthesis and measurement 
techniques^ have allowed to prepare nanostructructured 
materials with different dimesionalities on the submicron 
range. This has open the field of nanomagnetism to a 
handful of new opportunies^ that exploit new magnetic, 
optical and electrical properties that emerge when re- 
ducing the size of the particles to the nanoscale^^, of 
interest in wide areas of science ranging from magnetic 
recording and quantum computing^ to Earth sciences^ 
and biomedicine^ d ! 11 ! 12 ! 13 . 

Due to their reduced dimensions, nanoparticles display 
peculiar magnetic and transport properties*^ that are not 
present in the bulk materials as a consequence of the in- 
terplay between intrinsic properties arising from finite- 
size effects and collective effects due to different kinds 
of interparticle interactions^. A direct consequence of 
the finite size of the particles is superparamagnetism, 
which is a drawback for magnetic recording applications 
because it causes thermal destabilization of the record- 
ing units. However, superparamagnetic (SP) response 
is desirable for most of biomedical applications. Another 
effect influencing the magnetic response of the nanoparti- 
cles is the reduction of the magnetic net moment as com- 
pared to bulk. This is due to the competition between 
the different magnetic ordering at the particle core and 
its surface, which has a higher degree of disorder due 
to the broken symmetry, roughness and different stoi- 
chiometry from the bulk material. Particle surfaces are 
usually exposed to environment and are, therefore, easily 
oxydized, resulting in core/shell structures that can be 



otherwise produced by controlled chemical synthesi s 16 ' 17 
in a variety of morphologies and compositions. Magnetic 
core/shell nanoparticles with functionalyzed shells and 
coatings are also necessary in biomedicine for applica- 
tions in targeted delivery and diagnostics^. 

An attractive composition results from the combina- 
tion of a ferromagnetic (FM) core surrounded by an anti- 
ferromagnetic (AFM) shell (usually an oxide) coupled by 
the exchange interaction at the interface between them. 
Interesting proximity effects result from the structural 
modification and competition of different magnetic or- 
derings at the FM/AFM interface. In particular, the 
exchange coupling at a FM/AFM interface may induce 
unidirectional anisotropy in the FM below the Neel tem- 
perature of the AFM, causing a shift in the hysteresis 
loop, a phenomenon known as exchange bias (EB). For 
EB to occur, the Curie temperature Tq of the FM has 
to be greater than Tjy and the system has to be cooled 
from a starting temperatute in between in the presence 
of an applied field -Hfc- Moreover, the anisotropy of 
the AFM has to be high enough so that its spins remain 
fixed during the hysteresis loop. Although the first ob- 
servations of this phenomenon, dating back five decades 
ag o 19 i 20 , were reported on oxidized nanoparticles, most of 
the subsequent studies have focused on layered FM/AFM 
structure s 21 ! 22 because of their application in advanced 
magnetic devices^i 2 ^. However, in recent years, the study 
of EB in nanoparticles and nanostructures has gained re- 
newed interest since it has been shown that control of 
the core/shell interactions or of the exchange coupling 
between the particle surface and the embedding matrix 
can be a way to beat the SP limi t 25 ' 26 . 

Both nanoparticles and layered systems display com- 
mon phenomenolgy although, in the later wider 
range of experimenal techniques have been used which 
have provided deeper knowledge on the microscopic 
mechanisms that are at the basis of the EB effect. Thus, 
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knowledge of the magnetic structure at the interface has 
become a subject of primary interest in understanding 
EB. At difference from layered systems, the interface 
of core/shell nanoparticles naturally incorporates rough- 
ness and non-compensation of the magnetization, two of 
the main ingredients for which different assumptions are 
adopted by the existing models for EB in films 27,28 . How- 
ever, the interpretation of the results may be hindered 
by collective effects and interactions with the embedding 
matrix since, up to date, no EB experiment has been con- 
ducted on a single particle, which would allow to confront 
the results with the existing models. 

Most of the theoretical framework for the explana- 
tion of EB is based in macroscopic or phenomenologi- 
cal models for layered systems, adapted to the partic- 
ular structure and composition of specific combinations 
of materials. Guided by simplicity and reproducibility 
of experimental results, simplifying assumptions about 
the magnetic order in the FM and AFM layers are often 
assumed which may not allow to understand the real mi- 
croscopic origin of the EB effects. Moreover, and despite 
the similarities in both cases, the models used for EB in 
layered systems are not well suited for particle systems, 
since surface effects and the reduced dimensionality of 
the nanoparticles are supposed to play a role in the ob- 
servation of EB. 

For this purpose, computer simulations based either 
on Monte Carlo (MC) methods or on the micromagnetic 
approach 2 ^ have proved useful to gain insight into the 
microscopic origin of EB. These methods allow to take 
as inputs microscopic parameters such as exchange and 
anisotropy constants specific to the materials at hand and 
also to take into account the specific arrangement of the 
magnetic atoms in a lattice. As an output, macroscop- 
ically measurable quantities, such as the magnetization, 
can be computed without loosing valuable information 
about the microsocopic magnetic configurations that are 
at the origin of the observed phenomena. 

In this article, we will review the main phenomenol- 
ogy associated to EB in core/shell nanoparticle systems 
and the main existing models to explain it. The review 
is organized as follows. In Sec. II, we review the main 
nanoparticle systems for which EB has been reported, 
with special emphasis in nanoparticles with core/shell 
structure. Next, in Sec. Ill, we present a summary of 
the different phenomenology associated to EB found ex- 
perimentally for core/shell nanoparticles. In Sec. IV 
the results of MC simulations of a model of core/shell 
nanoparticle recently proposed by u o 30 ' 31 , together with 
other results in the literature, are presented. We end 
with the final conclusions and remarks in Sec. V. 



II. CORE/SHELL NANOPARTICLES 
DISPLAYING EB 

Observation of EB in nanoparticles has been reported 
for a wide variety of materials and morphologies which 



can be divided in three categories: (1) single phase ferri- 
magnetic or antiferromagnetic oxides, (2) nanoparticles 
embedded in a AFM matrices and (3) nanoparticles with 
core/shell structure. 

In the first group we have ferrites, manganites and an- 
tiferromagnetic particles (see tables 1 and 2 for a sum- 
mary of results in Ref. 124 ). The origin of EB in this 
kind of nanoparticles is not established yet, since, in this 
case, one cannot strictly speak of a FM coupled to an 
AFM material. Therefore, the observation of EB has 
been attributed to the freezing of a spin-glass surface 
layer of spins which is formed due to finite-size and sur- 
face effect o 32 ' 33 . However, the fact that for this kind of 
particles high field irreveribilities and non-saturating hys- 
teresis loops are commonly found, poses the question of 
whether minor loop effects could also be at the origin of 
loop shifts. 

Another way to get a high density of interface cou- 
pling FM and AFM phases of different materials is by 
embedding FM particles in AFM hosts synthesized with 
different techniques, although in these systems no clear 
separation between core and a well defined shell can be 
made. We refer the reader to Refs. (section 3.3) and 
[34] for recent reviews of results in this kind of systems. 

Finally, some of the largest observed EB fields have 
been reported for particles consisting of a FM core 
and an AFM (or ferrimagnetic) shell which has been 
grown around the core by chemical modification (usu- 
ally partial oxidation) of the FM material. Among 
them, some particular combinations have deserved 
special attention as C0/C0O, where the EB effect 
was first described by Meiklejohn and Bea n 19 ' 20 ' 35 
and revisited some decades later by Gangopadhyay et 
al i 36 i 37 and later by Peng and coworker a 38 i 39 i 40 i 41 i 42 . 
More recent studies of EB phenomenology in C0/C0O 
nanoparticles are listed in what follows [25! j 2 6lJ 4 3l j 44j . 

Other core/shell particle systems having Co and other 
oxides have also been studied such as NiCo/NiCoO 63 , 
Co/CoN 64 , Co/MnO^, Co 80 Ni 2 o/oxide^>^ and 
CoPt/CoO^ 8 -. Studies of iron oxidized particles such 
as Fe 3 4 /FeO^, Fe/j^Oj^^^^e^rs^ 

and Fe/FezO^^ 8 ^" 8 ^ 4 -" 8 ^ have also reported a 
variety of effects related to EB. Let us mention also 
the cases of Ni/NiO^^^^^^^ , Cr 2 3 /CrO:A 
ZnFe 2 4 /CoFe 2 04^, FeCo/CoFeO^., and also re- 
cently FePt/MnO^ 7 -, and FePt/FeaO^^. There 
has been also recent reports of EB in unconventional 
morphologies such as AFM Mn (core)/ferrimagnetic 
(FIM) Mn 3 4 (shell) 1 ^., FIM CoFe 2 4 (core) /AFM Mn 
(shell) 1 ^ 1 - nanoparticles, Fe304/Co nanocable o 102 ' 103 , 
and even Fe/Co oxidized particles encapsulated in a 
ferritin case 1 ^ 4 -. 
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III. EB PHENOMENOLOGY 

Although the main indication of the existence of ex- 
change bias is the observation of shifted hysteresis loops 
along the field axis after field cooling across the Neel tem- 
perature of the AFM T/v , some other macroscopic effects 
usually accompany the observation of loop shifts. In what 
follows, we will summarize the main experimental obser- 
vations related to EB peculiar to core/shell nanoparticles, 
comparing them with similar results observed in layered 
systems when possible. 



A. Coercivity increase 

The most usual is the increase in the coercive field He 
after field cooling observed below T/v, which is related 
to the unidirectional anisotropy induced on the FM by 
the field cooling process. Increased coercivities should 
appear only when the anisotropy of the AFM component 
is small compared with the exchange coupling with the 
FM component. In this case, partial rotation of the spins 
of the AFM shell, which are dragged by the FM core 
spins during the hysteresis loop, is expected resulting in 
increased He- A two times increase was found in for 13 
nm Co/CoO in RefM and also in Ref.- 1 ^. 



B. Partice size dependence 

As in the case of thin film systems, where the exchange 
bias field H e b is found to depend both on the thickness 
of the FM and AFM layers, EB effects in nanoparticles 
should depend on the particle size (core diameter Dq) 
and the thickness of the AFM shell Dgh- The dependence 
of H e b on the particle core size should be similar to that 
on the thickness of the FM layer in thin film systems and, 
therefore, H e b should increase when reducing the particle 
size H e b ~ 1/Dc- This trend was first reported in oxide 
passivated Co particle o 36 ' 37 in the size range of 5-35 nm 
and later confirmed by Peng and coworker a 38 i 39 i 41 i 42 on 
oxide coated Co/CoO particles with sizes 6-13 nm obtain- 
ing EB fields as large as 10.2 kOe for the smallest particles 
and a coercivity of 5 kOe. This has also been observed 
in oxygen passivated Fe particles with diameters of 6-15 
n m 71 ' 72 , and in Fe/7-Fe203 particle a 76 ' 79 . Moreover, a 
critical particle size below which EB is absent for any ra- 
tio of ferromagnetic and antiferromagnetic constituents 
has been reported^ for Co/CoO 3 nm nanoparticles em- 
bedded in AI2O3. The reason is that, due to the large 
surface-to-volume ratio below the critical size, the ex- 
change energy at the FM-AFM interface becomes smaller 
than both the effective Zeeman energy of the FM and the 
anisotropy energy of the AFM. In a later study on a sam- 
ple consisting of 2.5 nm Co clusters embedded in a CoO 
matrix, the same authors performed a more complete 
study for samples with different oxide layer thicknesses^ 7 -. 



In fact, some authors have also reported an upper criti- 
cal size (40 nm for the CoNi/CoO partices embedded in 
PVC of Ref. [13) for the observation of EB. In another 
study of CoFe204 particles 1 - ^ with diameters 15-48 nm 
a nonmonotonic size dependence of H e b , similar to what 
is observed in He , has been observed (although at much 
higher temperature of 77 K), with an increase with par- 
ticle size up to a peak at around 27 nm followed by a 
subsequent decrease and vanishing for 40 nm particles. 
More recently, a study by Boubeta et aJ.— on oxidized 
Fe particles with diameters ranging from 5 to 13 nm have 
confirmed the disappearence of EB below a critical diam- 
eter of 5 nm and attributed this effect to the decreasing 
thickness of the spin-glass-like layer when decreasing the 
nanoparticle size. However, in the oxidized Fe particles 
studied by Ceylan et aii£&, the small particles (7.5 nm 
in diameter) were found to have much higher Heb than 
the big ones (13 nm in diameter), probably due to the 
increased relative effect of the AFM shell and the more 
amorphous structure of the shell in the smallest particles. 



C. Shell thickness dependence 

Fewer studies have focused on the role played by 
the shell thickness, since the formation of oxidized 
phases cannot be easily controlled independently of the 
core size. As indicated by some models of EB for 
thin ^maiSS^SMmUMUmiUli, there should be 
a minimum critical shell thickness for the observation of 
a loop shift, since the anisotropy energy per unit area of 
the AFM has to be larger than the interfacial exchange 
energy for EB to exist. Above this limiting thickness, 
H e b should increase with D$h up to a critical shell thick- 
ness above which it would become independent of Dsh- 
This has been partially corroborated by several works 
on nanoparticles of different compositions in which sam- 
ples prepared by the same technique but different degrees 
of oxidation were compare d 51 ' 55 ^ 7 -^ 3 -^ 1 -'^^ 3 -^ 4 -'^. More- 
over, the critical shell thickness in nanoparticles should 
depend on the anisotropy of the AFM as was first estab- 
lished in bilayers by Lund et alf 1 ^. 



D. Training effects 

A less studied effect, first described for thin films 1 -^ 7 - 
but also observed in nanoparticles, is the so-called train- 
ing effect, which is observed when the hysteresis loop are 
succesively repeated a number of times n after FC. H e b 
gradually decreases with n in thin films, reflecting the de- 
viation of the AFM spin structure at the interface layer 
from its equilibrium configuratio n 118 ' 119 . The relaxation 
of the frozen spins along the cooling field direction re- 
duces the effective pinning energy, resulting in a decrease 
ofH eb with the number of field cycles. Moreover, the bias 
field increase with increasing sweep rate of the magnetic 
field has been described by a dynamically generalized the- 
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ory based on triggered relaxation, in excellent agreement 
with the experiment o 120 ' 121 ' 122 . A quantitative explana- 
tion based on the Kolmogorov-Avrami model describing 
the dynamics of AFM layers 1 ^ 3 , seems to describe cor- 
rectly experimental data on the H e b(t) dependence. Also 
the symmetry of the anisotropy in the AFM seems to be 
crucial for the understanding of training effects 1 ^. 

In core/ shell nanoparticles, this training effect is char- 
acterized by a decrease of the coercive field on the de- 
scending field branch of the loop, whereas the ascending 
branch is usually retraced on succesive cyclings. More- 
over, the training rate seems to depend strongly on the 
properties, namely the AFM or ferrimagnetic character, 
of the oxide shells. Thus, whereas in some Co/CoO parti- 
cle systems the training is more pronounced after the sec- 
ond cycl o 41 ' 57 ' 61 , in some Fe/Fe oxide particle systems 8 - 2 -, 
the training effect is only decreas ed to about 89% after 
the 14 th cycle (see also Refs. I76lfl25h. Clearly related 
to training effects is the observation of aging effects on 
H e b when the hysteresis loops are measured at different 
waiting times after the cooling field is applied^. 



E. Temperature dependence 

Of course, both H e b and H c are thermal dependent 
quantities. Since the AFM or ferrimagnetic magnetic or- 
der at the particle shell, which is at the origin of the 
existence of EB, is degraded by temperature, EB should 
disappear when approaching the ordering temperature of 
the shell T/v, which is lower than the Curie temperature 
Tc of the FM core. In fact, for most experimental sys- 
tems, EB disappears at a so-called blocking temperature 
Tb lower than T/v, although this is not necessarily true 
for H c , for which finite values higher than those obtained 
after ZFC are usually observed up to Tv 1 ^ 7 -. 

For thin films, it has been argued that the difference 
between Tb and Tv depends on the AFM layer thickness 
and is not related to finite-size effects on T/\42i. However, 
in particle systems, this has been attributed to the SP 
behavior of the AFM oxide shell at a temperature lower 
than the T/v of the shell, which might be composed of 
very small crystalrite a 37 ' 45 ' 48 . With respect to the exact 
T dependence, in thin films, linear dependencies of both 
quantities are usually observe d 127 ' 129 ' 130 in accordance 
with the random field model of Malozemof f 131 ' 132 ' 133 . 
This is not always the case for core/shell nanoparti- 
cles, for which faster than linear decays of iTb and 
H c have been reported for Co/CoO particle s 37 ' 41 ' 48 , al- 
though quasi-linear dependencies are also found^ 9 -. A 
law of the kind H eb (T) = H eb (Q)(l - T/T N ) n with 
n = 3/2 has been show n 43 ' 45 ' 47 to fit experimental data 
on Co/CoO particles, which is in accordance with the 
predictions of a model for polycrystalline bilayers 1 ^ 4 - that 
takes into account the thermal instability of the AFM 
shell. It must also be remembered that, when dealing 
with nanoparticle systems, other factors apart form the 
structural ones, intrinsic to the particle, such as the vol- 



ume distribution, randomness of the anisotropy axes and 
the existence interparticle interactions 1 * may influence the 
thermal dependence of both H e b and H c . 



F. Cooling field dependence 

There is no general trend for the dependence of H e b 
on the cooling field magnitude in layered systems. De- 
pending on the details of the microscopic structure of the 
interface and the AFM layer and the preparation con- 
ditions, both a slight decrease^ or increase^ of H e b 
with increasing T have been reported. However, some 
systems 21 display loops shifts towards positive field val- 
ues instead of to negative fields for large cooling fields. 
This effect has been argued to be possible when the 
coupling at the interface is AFM. Estimations of the 
crossover field have been give n 136 ' 137 ' 138 ' 139 , and experi- 
ments have also proved the validity of the hypothesis in 
several bilayered systems^^Iii^i 2 .^^. 

Field cooling dependencies have been reported only 
recently in core/shell nanoparticles. In CoFe2 04 
particles 1 ^, H e b has been found to increase with the cool- 
ing field for values of Hpc up to 5000 Oe while, for higher 
fields, a slight decrease is observed accompanied by a de- 
crease in the vertical loop shift. On the other hand, while 
for Co/CoO nanoparticles^ 9 - H e t continue to increase for 
fields up to 5 T with values of the order of 1-2 kOe at 
300 K, for Fe/FeO nanoparticle o 70 ' 73 ' 145 , H eb presents a 
maximum at a field cooling value around 5 kOe which 
increases with decreasing T. In this case, H e b decreases 
with further increasing the cooling field, reaching a value 
of only 250 Oe at 5 T and 5 K. The authors argued 
that the appearance of the maximum is due to the glassy 
magnetic nature of the oxide phase at the shell, which 
might be destroyed by increasing magnetic fields or tem- 
peratures. A similar behaviour has also been reported in 
phase-separated LSCO perovskiteA 4 - 6 . A clear-cut inter- 
pretation for these systems is still lacking. 

G. Asymmetry of the hysteresis loop 

Another commonly observed feature in bilayers is an 
asymmetry between the descending and increasing field 
branches of the loops after FC, which has been related 
to different magnetization reversal mechanisms in each 
of the branches. While in the descending field branch 
reversal takes place usually by uniform rotation, in the 
increasing field branch, reversal by nucleation and prop- 
agation of domain walls or non-uniform structures seems 
to be the dominant mechanism. Different techniques, 
sensitive to microscopic magnetic configurations of the 
FM and AFM, have confirmed these different reversal 
mechanisms. First studies on this issue were performed 
in FeNi/FeMn films by magneto-optical method o 147 ' 148 
and in MnF 2 /Fe 1 i 9 -^ and CoO/Co bilayer o 127 ' 151 ' 1 ^ 2 .' 1 ^ 3 . 
by polarized neutron reflectometry. Later on, also X-ray 
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photoemission microscopy has been used in Fe/MnPd 
films 1 ^ 4 -, time-resolved Kerr magnetometry in FeF2/Fe 
bilayergA5fi and neutron scattering in patterned Co/CoO 
nanostructures 1 ^. Recently, the origin of asymmetric 
loops in some particular systems has been ascribed to 
the competition between the FM and the interfacial FM- 
AFM exchange anisotropie o 157 ' 158 and different reversal 
processes in both loop branches have also been revealed 
by SXRMS and techniques in perpendicularly coupled 
exchange coupled films 1 ^ and by MOKE in Fe/MnF 2 
bilayers 1 ^ -. 

Asymmetries in the shape of the hysteresis loops of 
core/shell particles are also evident in some systems, but, 
in this case, present experimental techniques cannot eas- 
ily give information about the microscopic mechanisms 
involved in the reversal processes because of the particle 
size dispersion always present in samples. For this pur- 
pose, experiments being able to measure magnetic prop- 
erties of a single nanoparticle (in the spirit of those per- 
formed by Wernsdorfer and coworker a 161 ' 162 ' 163 ) would 
help to clarify this controversial issue. 



H. Vertical loop shifts 

In some systems, shifts along the magnetization axes 
have also been reporte d 140 ' 164 that have been related 
to induced magnetic moments. This vertical shift de- 
pends on the cooling field (it may be negative for low 
Hfc and positive for large Hfc) and the microstruc- 
ture of AFM layer. Recently, X-ray magnetic circu- 
lar dichroism (XMCD) experiments on Ni/FeF2 bilay- 
ers have proved that the vertical shift is due to the 
existence of uncompensated Fe pinned moments in the 
AFMl 4 ^. Huang and co-workers 1 ^ 4 - have observed lin- 
ear dependence of the exchange field on the magnetiza- 
tion shift in ZnCoo-070/NiO layers, proving the role of 
uncompensated pinned spins on the observation of the 
effect. Some core/shell nanoparticles also display this 
phenomenology. Vertical shifts have been reported*^ for 
Ni/Ni02i>M, Co/CoO^se^^ 1 -, Fe/Fe 2 0^^- and in 
milled Fe/MnOz 1 ^ particles, with values much higher 
than those reported for bilayers. The linear dependence 
of the vertical shifts measured at different temperatures 
on H e f, found in I49II75I indicates that the vertical shifts 
are proportional to the number of net frozen spins. A 
nonmonotonic dependence of the shifts on the particle 
size and cooling field, in agreement with that found for 
Heb, has been reported by Mumtaz et al.^2&. 

I. Nature of the interface coupling 

Recently, several spectroscopic techniques have pro- 
vided insight on the structure and magnetic behavior 
of the interface spins at a microscopic level, demon- 
strating the crucial role played by uncompensated in- 
terfacial spins on EB in several bilayered thin film 



system a 154 i 159 i 166 i 167 i 168 i 169 i 170 and also demonstrating 
unambiguously the existence of domain walls in the FM 
parallel to the AF/FM interface 1 - 7 -!'. Similar techniques 
applicable to nanoparticles such as X-ray absorption and 
XMCD have also been used to study Fe oxide passivated 
iron nanoparticles 1 ^. The relative sign of the metal and 
oxide related dichroism allows to conclude that the cou- 
pling across the interface is FM. This finding is opposed 
to the situation at the Fe(110)/Fe3O4 interface, where 
an AFM coupling was found 17 - 3 -. Presence of uncompen- 
sated Co magnetic moments at the interface of a 2U2.5 
nm CoO shell surrounding a metallic fcc-like 7-8 nm Co 
core was also evidenced by XMCD 44 ' 47 . 



J. Other recent observations 

In this last subsection, we would like to mention 
some very recent experimental observations in core/shell 
nanoparticle systems which have given evidences of new 
phenomenology not mentioned in the previous subsec- 
tions and that we think will estimulate further stud- 
ies both from the experimental and theoretical point of 
views. Tracy et al£L have reported an investigation of 
the role of defects on the magnetic properties of Co/CoO 
nanoparticles in which, by measuring magnetization and 
thermoremanence curves under ingenious FC protocols 
with intermediate field reversals, they are able to show 
that the defect moments freeze at low temperature and 
have a distribution of melting temperatures and that they 
dominate EB at low temperature, exhibiting also a ther- 
mal memory effect. The role of dilution on the AFM 
have also been studied in bilayer a 174 ' 175 . Both experi- 
mental and simulation results confirm an enhancement 
of H e b with increasing defect concentration. Nogues et 
al— have demonstrated that the magnetic properties of 
Co/CoO nanoparticles embedded in an AI2O3 matrix, 
depend strongly on the in-plane coverage, even in the di- 
luted regime. In particular, the authors have found that 
both He and H e b radically increase with increasing cov- 
erage. The experiments allow the authors to conclude 
that these observations cannot be accounted by dipolar 
interactions between the cores and should be attributed 
to shell mediated interactions when particle become in 
contact. This would also help to explain the scatter of 
values for He and H e b found in the literature. A study of 
CoO granular films deposited on layered FM structures 
by Gruyters^ 2 - have shown that EB in this system can be 
explained by the spin-glass-like state in the nanoparticles 
constituting the CoO film without the need for core/shell 
structure. These results show that pinning effects in EB 
systems are not only related to uncompensated spins, 
but may arise due to a frozen state in the AFM similar 
to a spin-glass. Moreover, the deduced unusually large 
uncompensated magnetization has no simple quantita- 
tive relation to H eb , a fact that requires further theoret- 
ical developement in order to be understood. The same 
author has proposed a modeU^ 6 ., based on the random 
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magnetic anisotropy of CoO nanoparticles, according to 
which the observation of EB can be attributed to an in- 
teraction between the AFM order and uncompensated 
spins in the AFM material without explicitly invoking 
the exchange coupling to a FM. Another study that will 
hopefully provide a new direction for studies of EB is that 
by Ali and co-workers^ 7 - on a Co/CuMn bilayer system, 
which has evidenced the possibility of observing most of 
the EB associated phenomenology using a spin-glass ma- 
terial instead of conventional AFM. A striking difference 
form FM/AFM bilayers a change in sign of the bias field 
just below the blocking temperature has been found in 
this system, indicating that the indirect RKKY exchange 
within the pinning layer may account for the observed ef- 
fects. One may wonder if core/shell particles with similar 
morphologies could also give surprising new effects. 



IV. MODELS AND SIMULATIONS 

Some microscopic models for bilayers have un- 
dertaken calculations of EB fields under certain 
assumption o 178 ' 179 , numerical studies based on 
a mean field approach^ or Monte Carlo (MC) 
simulation a 181 i 182 i 183 making different assumptions 
about the interface. However, only very recently, some 
works partially addressing the EB phenomenology in 
nanostructures have been publishe d 184 ' 185 . 



A. Model of core/shell particle 

In order to understand what is the microscopic origin of 
all the phenomenology associated to EB effects presented 
in the preceding section, we have developed a model for a 
single nanoparticle with core/shell structure which cap- 
tures the main ingredients that are believed to be neces- 
sary for the observation of EB. A schematic drawing of 
the particle is shown in Fig. [TJ Atomic spins are consid- 
ered to sit on the nodes of a sc lattice and the particle is 
buildt by considering the spins inside a sphere of radius 
R (measured in multiples of the unit cell dimensions a) 
centered in on of the lattice nodes. Three regions are 
distinguished inside the particle: a core with radius Rc, 
a shell of thickness Rsh = R — Rc and the core/shell 
interface that is formed by the core (shell) spins hav- 
ing nearest neighbours on the shell (core). In most of 
the results presented in the following, we have consid- 
ered a fixed particle size R = 12a an a shell of thickness 
i?Sh = 3a. Taking a — 0.3 nm, such a particle corre- 
sponds to typical real dimensions R ~ 4 nm and i?sh — 1 
nm and contains 5575 spins, of which 45 % are on the 
surface. Since we are interested in studying magnetic 
properties observed in real core/shell particles, we will 
consider that the core of the particle is made of a FM 
material and that the outer shell is an AFM. Different 
characteristic microscopic parameters, such as exchange 
and anisotropy, will be considered in the three regions, 




FIG. 1: (Color online) Schematic drawing of model of a 
core/shell nanoparticle of total radius R used in the MC sim- 
ulations. The spins sit on the nodes of a sc lattice. The AFM 
shell has width Rsh (green and yellow spins) and the FM core 
(blue spins) a radius Rc = R — Rsh- The core/shell inter- 
face (light blue and yellow spins) is formed by the core (shell) 
spins having nearest neighbours on the shell (core). 



with fixed values at the core and shell regions and that 
will be varied at the interface in order to study what is 
its specific role in establishing EB properties. 

To account for the finite values of anisotropy in real 
systems, we have considered a model of Heisenberg clas- 
sical spins Si, interacting according to the following mi- 
croscopic Hamiltonian 



H/ki 



-Jc ^ S t ■ Sj - Js 22 * ' 

(ij)ec (i,j)esh 

-Jmt S t ■ Sj ~ k c Y,(SD 2 

(<ec,j'eSh) »ec 

N 

ieSh i=l 



(1) 



The first three terms describe the nearest-neigbor ex- 
chage interactions between the spins with different values 
of the exchange constants at the different particle regions. 
Core spin are FM with Jc > 0, whereas spins in the shell 
are AFM with Js < 0. The values of these constants will 
be kept constant and fixed arbitrarily to Jc = 10 and 
Jsh = — 0.5Jc, which just fix the Curie temperature of 
the FM to T c = 29 K and the Neel temperature of AFM 
to Tn — 14.5 K, a value lower than Tc as is the case in 
most oxides with respect to their native materials. Fi- 
nally, for spins the exchange constant at the interaface 
Jint ^ will be allowed to vary between and ±Jc in 
order to study the role played by the coupling across the 
core/shell interface on magnetic properties. 

The fourth and fifth terms correspond to the on-site 
uniaxial anisotropy with k c and k$h the values of the 
anisotropy constants at the core and at the shell. They 
can be related to values in real units through the corre- 
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T(K) 



FIG. 2: (Color online) Thermal dependence of the normalized 
magnetizations of a core/shell particle when cooling form a 
disordered state at T > Tjv down to T = 0.1 in the presence 
of an external magnetic field h?c = 4 K. The values of the 
exchange coupling at the interface are (a) Jint = — 0.5Jc and 
(b) Jint = +0.5Jc. The different curves correspond to the 
contributions of the core Mc, shell Msh and interface Mint 
spins to the total magnetization Mr. Insets display the con- 
tributions of only the interfacial shell spins to Mint (Mgf , in 
circles) and, among these, the ones having 1 (black), 2 (red) 
and 3 (green) nearest-neigbors in the core. 
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FIG. 3: (Color online) (a) Spin configuration of an equatorial 
cut of the particle attained after the field cooling process de- 
scribed in Fig. 1. Core spins are dark blue, spins at the shell 
are green while inteface core and shell spins have been col- 
ored in light blue and yellow, (b-d) Configurations of the core 
(shell) spins at the interface having 1, 2 or 3 nearest neigh- 
bors in the shell (core). [Reprinted with permission from Ref. 
EH, O. Iglesias et al. Physica B 372, 247 (2006). Copyright 
©Elsevier B. V.l 



Carlo simulations using Metropolis algorithm. As for the 
spin updates, we use a combination of the trial steps 
which has proved useful for Heisenberg spins with finite 
anisotropics as described elsewher o 188 ' 189 . 



spondence 



B. Field cooled states 



K C V K Sh S 



(2) 



, where Kq and K$h are the anisotropy constants in 
units of energy per unit volume (V) or surface (S) of 
the particle. The value of kc will be fixed to kc = 1 
K, which just sets the value of the anisotropy field of 
the FM core, whereas the anisotropy at the AFM shell 
has to be higher than that in the core as required to 
pin the AFM spins during the hysteresis loops so that 
EB is observed. Therefore, we fix fcsh = 10 K, which is 
which also in agreement with the reported enhanced sur- 
face anisotropies due reduced local coordination at the 
outer particle sheil a 186 ' 187 . Finally, the last term is the 
Zeeman energy coupling to an external magnetic field H, 
where h — fiH/ks denotes the field strength in temper- 
ature units, with \x the magnetic moment of the spin. 

Based on this Hamiltonian, we have performed Monte 



First, we will study the magnetic state of the particle 
after a field cooling process with the purpose to charac- 
terize the magnetic order induced on the interfacial spins. 
Our protocol to simulate the field cooling process is as 
follows. We start the simulations from a high temper- 
ature To > Tjv disordered state in which the spins are 
pointing in random directions with zero net magnetiza- 
tion. The temperature is then reduced in constant steps 
S = 0.1K down to the final temperature T = 0.1 K in 
the presence of a magnetic field five = 4 K applied along 
the easy-axis direction. At each temperature, the mag- 
netization is averaged over a number of 10000 MC steps 
after 10000 MC steps used for thermalization, using the 
usual heat bath dynamics for continuous spins. 

As an example, the thermal dependence of the nor- 
malized magnetization along the field direction is shown 
in Fig. [2] for a particle with AFM or FM interface cou- 
pling Ji n t = =f0.5Jc. In this figure, the contributions 
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h (K) h (K) 

FIG. 4: (Color online) Hysteresis loops for a particle with ra- 
dius R = 12 a obtained from a ZFC state and after FC down to 
T = 0.1 in a field /ifc for Jsh = — 0.5Jc and Jint = — (+)0.5 Jc 
in the left (right) column. Panels (a) display the total normal- 
ized magnetization component along the field direction. Pan- 
els (b) show the normalized contributions of the shell spins 
at core/shell interface to the total magnetization of the loop. 
Panels (c) show the contribution of the interfacial spins at the 
shell to Mjjjjf having 1, 2 or 3 nn in the core. 



of the spins in the core (Mc), in the shell (Msb.) and 
at the interface (Mi nt ) to the total magnetization Mr 
have been plotted separately. As it can be seen in the 
main panels, during the cooling process, the core spins 
progressively order ferromagnetically as indicated by the 
increase of Mq towards 1. At the same time, as T is re- 
duced below the Neel temperature of the shell, the AFM 
order is also established in the shell spins, although a fi- 
nite value of Msh remains at the lowest temperature due 
to the noncompensation between sublattices caused by 
the finite-size and spherical shape of the particle. Most 
importantly, independently of the nature of the coupling 
between the core and shell spins, the interfacial spins 
are not compensated, as indicated by the finite magne- 
tization attained at low T, which, of course, is lower 
in the AFM case (M Int = 0.37) than in the FM one 
(Mi nt = 0.605. 

In order to gain deeper understanding on the origin 
of this net interface magnetization, first notice that the 
interfacial spins at the core are all pointing in the field di- 
rection after the FC process. Therefore, uncompensated 
moments must be originated at the shell interface. We 
show in the insets of Fig. [2] the contributions to M Int 
of the interfacial spins at the particle shell in the curve 
labeled Mg£* (in circles) . Comparing the insets in panels 



(a) and (b) , we see that the sign of the net magnetization 
at the shell interface is in accordance with sign of the in- 
terface coupling. Further inspection of the contributions 
of spins having different number of nn in the core, allows 
us to conclude that the aligning effect of the cooling field 
is more effective for the spins with lower number of nn in 
the AFM case and for those with more nn in the FM case. 
It can also be seen that the major contribution to the net 
interface magnetization comes from the shell spins with 
3 nn in the core, independently of the sign of Ji nt . As 
it can be concluded from the preceding observations, the 
geometric structure and magnetic ordering of the inter- 
face in a core / shell nanoparticle is more intricate than in 
the case of FM/AF coupled bilayers due to the roughness 
inherent to the geometry of the interface. At difference 
from bilayers, interfacial spins may have different number 
of neighbors depending on their position and, therefore, 
the interface spins at the shell present regions with either 
local compensated or uncompensated magnetic order. 



C. FC hysteresis loops 

In order to study the phenomenology associated to EB 
effects, we have also performed simulations of hysteresis 
loops following a protocol that mimics the experimental 
one: configurations obtained at the lowest temperature 
after the FC process described in the preceding section 
are used as the starting state and then the hysteresis 
loops are recorded by cycling the magnetic field from h = 
4 K to h — —4 K in steps 5h = —0.1 K and the different 
quantites averaged during 200 MC steps per spin at every 
field after other initial 200 MC stepd are discarded for 
thermalization. Hysteresis loops obtained from a zero 
field cooled (ZFC) state have also been simulated starting 
from a demagnetized state at the measuring temperature, 
then following the first magnetization curve up to h = 4 
K and, finally, performing the hysteresis loop as described 
before. 

Typical ZFC and FC hysteresis loops are shown in Fig. 
H] (upper panels) for two values of the interface coupling 
Jint/Jc = —0.5, +0.5. Compared to the loops obtained 
from ZFC state, the loops obtained after FC are shifted 
towards negative field values and have slightly increased 
coercivity (see Figs. UK), independenly of the sign of the 
interfacial exchange coupling. The values of the coercive 
fields for the decreasing and increasing field branches will 
be denoted by and /ij, respectively. Therefore, the 
coercive field and the EB fields are defined as he = {h^ — 
/iq)/2 and h e b = (h^ + h^)/2, respectively. The origin 
of the shift in the FC case can be better understood by 
looking at the contribution of interface spins belonging to 
the shell, Mg^*, to the total magnetization as displayed 
in the middle panels of Fig. [H 

As we have previously revealed by the detailed in- 
spection of the microscopic configurations attained after 
FC, the interfacial spins at the shell acquire a negative 
(Jint < 0) or positive (Ji n t > 0) net magnetization after 
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FC, in both cases higher than the one attained after ZFC, 
although more pronounced for the negative coupling case. 
This net magnetic moment, induced by the geometrical 
symmetry breaking and the alignment of groups of spins 
into the field direction, generates local fields on the core 
spins that point into the same direction as the external 
field, causing the shift of the hysteresis loops. 

To further support this observation, we note that the 
hysteresis loops are shifted by the same amount but to- 
wards the positive field axis when cooling in a field ap- 
plied in a direction negative with respect to the measur- 
ing field (see for example the dashed lines in Fig [6] for 
Jint/Jc = —0.5,-1). These observations reflect that, 
after the FC process, a fraction of the interfacial spins 
(«15 % of the interface spins at the shell) have been 
pinned along a direction compatible with the core/shell 
exchange interaction, as corroborated also by the verti- 
cal shifts in the Mg£* loops (to be commented below). 
This is no longer true for the ZFC case, for which a high 
fraction of interfacial spins follows the reversal of the FM 
core, as reflected by the change in sign of M§£ along the 
hysteresis loop. Moreover, FC hysteresis loops obtained 
for the same particles but without increased anisotropy 
at the AFM shell (performed setting ksh = 1) display no 
EB but, instead, have increased coercive fields compared 
to ZFC loops. In this case, no interfacial shell spins are 
pinned and, during reversal, they are dragged by the core 
spins due to the dominance of exchange coupling over 
anisotropy energy. This observation demonstrates that 
high anisotropy AFM are required to obtain exchanged 
biased loops. 

It turns out that disorder and frustration at the sur- 
face induced by radial anisotropy and finite-size effects 
alone are not enough to produce sizable loops shifts 
as simulations performed for particles with no AF shell 
demonstrate^. 



D. Quantifying h^: microscopic origin of EB 

One of the most controversial points in EB concerns 
the evaluation of the loop shifts from a model of the 
system at hand. Different theories and models usually 
predict EB shifts that differ by orders of magnitude from 
that measured experimentally. An archetypical example 
is the expression first derived by Meiklejohn and Bean 
(MB)±&2&2£ for a bilayer that reads 



Jeb 



/j. M F t F 



(3) 



where J e b is the interfacial exchange energy per unit area 
and Mf, tF are the magnetization and thickness of FM 
layer, respectively. Although this expression describes 
correctly the linear decrease of H e b with tp, it fails in 
the quantitative prediction of most of the measured loop 
shifts, the reasons being, essentially, that the FM/AFM 
interface is supposed to be fully uncompensated and ide- 
ally smooth and that the AFM is considered to be single 



domain with spins that remain unchanged during the re- 
versal of the FM. Other models based on refined versions 
of the MB model gave improved expressions for H e f, that 
agreed more reasonably with experimental values in some 
layered systems. Let us briefly recall that the model by 
Malozemo ff 131 ! 132 ' 133 , that incorporated the roughness of 
the interface as a random field acting on the FM layer 
and a model by Maurii^I that, following the pioneering 
work by NeeU^ 2 -, accounted for the possibility of domain 
wall formation in the AFM, gave modified expressions 
for the EB field of the kind H eb ~ M f//^ F ( A(J af being 
the domain wall energy density in the AFM), which re- 
sult in reduced values with respect to the MB model (see 
also the models by Kiwi at a h 28 i 138 i 139 and Stamps and 
co-worker a 27 i 193 ) . 

In spite of the profusion of models presented above, 
none of them takes into account the evolution of the spin 
structure of the FM and the AFM along the hystere- 
sis loops and this is the reason for their lack of agree- 
ment with experiments. More microscopic approaches 
such as the work by Takano et ali^, in which, by cal- 
culating the density of interfacial uncompensated spins 
in permalloy/ CoO bilayers, the authors predicted the 
correct magnitude of the exchange field, as well as the 
observed inverse dependence on interfacial grain size, 
have been more successful. More recently, a semi- 
quantitative account of the EB field magnitude has been 
presented in a simplified model for Co nanoparticles em- 
bedded in a CoO matrix-^. In order to link the mea- 
sured loop shifts to the microscopic details of the sam- 
ples, Monte Carlo and micromagnetic simulations based 
on microscopic models 19 - have proved useful. Among 
them, let us mention here that, to our knowledge, only 
the domain state (DS) model proposed by Nowak and 
collaborator a 130 i 181 i 196 i 197 i 198 have been able to establish 
a numerical correspondence between H e b and microscopic 
parameters by proving that H e b is proportional to the 
irreversible domain state magnetization of the AFM in- 
terface layer mi ds as if e {, = Jlnt,n ' DS where I is the FM 
layer thickness and \x the atomic magnetic moment. 

In the case of a core/shell nanoparticle, the analysis is 
more intricate due to the peculiarities of the core/shell 
interface as already commented in Sec. IIV B} and a more 
detailed analysis is needed. In order to elucidate the 
role played by the interface in establishing the EB ef- 
fect, we have studied the variation of h^,, /ij, he and /i c b 
with the interface exchange coupling J\ n t, presented in 
Fig. EH for negative Ji nt values. With increasing Ji nt , 
both tig and /ij decrease in absolute value, although 
they seem to reach a constant value when approaching 
I Jint | = Jc- As a consequence, a decrease in he and 
an increase in h e b is observed, with a nearly linear de- 
pendence, at least for values of |Ji n t| smaller than the 
exchange coupling at the shell Jgh = — 0.5 Jo Similar 
linear dependencies have been found in the DS model 
and some other models of bilayersA2&. MC simulations of 
a cylindrical nanodoli^, also demonstrated an increase 
in h e h with the scaled effective unidirectional anisotropy. 
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FIG. 5: (Color online) (a) Variation of the coercive fields h^, 
/iq, /ic and the exchange bias field h e b with the exchange 
coupling constant at the core/shell interface for Ji, lt < for a 
particle with Jsh = — 0.5Jc. (b) Variation of /i e j, with JW < 
(open circles) and JW > (open squares). The exchange bias 
fields computed from Eq. \E\ as described in the text are shown 
as filled symbols for Ji n t < (down triangles) and Ji n t > 
(up triangles). 



For both Ji n t ^ 0, the values of he and /i e b are very 
similar, as can be seen in Fig. Eb- With the increase of 
| Jint |, core spins become more coupled to the unpinned 
shell spins, therefore facilitating the magnetazation re- 
versal with the subsequent decrease in the coercivity, an 
observation also found in micromagnetic simulations of a 
model of coupled bilayers with grains in the AFM, which 
exhibit random uniaxial anisotropy and are weakly ex- 
change couple d 200 ' 201 ' 202 . At the same time, increasing 
| Jint | while keeping the values of Jc, Jsh and hpc con- 
stant, results in higher local exchange fields created by 
the uncompensated spins at the interface, causing an in- 
crease of the loop shift. Notice, however, that, increas- 
ing | Jint | above Jsh do not result in a further increase 
of h c \y, which seems to converge to a common value for 
both cases. The reason for this departure from linearity 
stems will be commented in the next sections. Finally, 
let us also mention that the values of the coercive and ex- 
change bias fields obtained from simulations are within 
the correct order of magnitude when expressed in real 
units. For example, for Ji n t/Jc £ [—0.3,-0.5], we ob- 
tain H c » 1.3 - 1 T and H ch w 0.27 - 0.43 T, which 
are in agreement with typical values found in studies of 
oxidized nanopartides 2 ^^ 3 -^^^! 7 ^. 

The proportionality of h ch to Ji nt should be taken as 
a hint for the microscopic origin of the loop shifts. As 



we have mentioned in previous paragraphs, the observed 
vertical displacements of the loops corresponding to the 
interface shell spins point to the existence of a net magne- 
tization at the core/shell interface due to uncompensated 
pinned spins at the shell interface^. If this is the case, 
the coercive fields after FC can be thought as the sum 
of the ZFC coercive field h^, and the local field acting on 
the core spins due to the net interface magnetization of 
the shell spins, so that they may be computed a a 30 ' 179 



h± = h° c 



JlntAfl 



Int. 



(4) 



where M± t = £ j£ 

{int,sh} z i$i is the net magnetization 
of the interfacial shell spins at the positive (negative) co- 
ercive fields , and zi is the number of nearest neighbors 
of spin i. Therefore, the coercive and exchange bias fields 
can be written as 



h c = h° c + J lnt (M+ t -M l - t )/2 
h eh = Jmt(M+ t +M- t )/2. 



(5) 



These expressions establish a connection between the 
coercive fields and loop shifts observed macroscopically 
and microscopic quantities that, although may not be 
directly measured in an experiment, can be computed 
independently from the simulation results. 

The values of /i c b obtained by inserting the M^ t val- 
ues extracted from the Fig. [4Jd in Eq. [5] are represented 
as filled symbols in Fig. [5b, where we can see that the 
agreement with the h c h values obtained from the hys- 
teresis loop shift is excellent within error bars. Recent 
experiments by Morel et ali 2 ^ 3 - on Co particles embed- 
ded in MnPt films have observed a clear correspondence 
between the measured H e b and Mafm, the magnetiza- 
tion induced in the AFM MnPt established by suitable 
FC procedures, which reinforce the validity of our model. 
Only for | Ji n t| > Jsh, an increase in | Ji n t| does not result 
in a further increase of /i c b, as reflected by a departure 
from linearity implied by Eq. [5j which means that the 
interfacial net magnetization at the shell may be acting 
on core magnetization components transverse to the field 
direction. 



E. Loop asymmetries 

In addition, a clear asymmetry between the upper and 
lower loop branches developes when increasing the value 
of the interface coupling, as it is apparent when compar- 
ing the descreasing and increasing branches of the loops 
in the top panels of Fig. [4^ or Fig. [6l This feature 
can be more clearly seen by looking at the average abso- 
lute value of the magnetization projection along the field 
axis through the reversal process, M% — Y^i \&i ' ^li dis- 
played in the middle panels of Fig. [6] for the core spins. 
This quantity presents peaks centered around the coer- 
cive fields that indicate deviations of the core magneti- 
zation from the applied field direction. In the ZFC case, 
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FIG. 6: (Color online) Upper panels display the hysteresis 
loops obtained after FC down in a field Tifc = 4 (circles) and 
Iifc = —4 (dashed lines) for three values of the exchange cou- 
pling constant Ji n t at the core/shell interface. Lower panels 
show the average magnetization projection of the core spins 
along the field axis m„ (squares) and the hysteresis loops for 
the component of the magnetization transverse to the field 
dire ction M tr (circles). [Reprinted with permission from Ref. 
l21ll . O. Iglesias et al. J. Magn. Magn. Mater, (in press, 
doi:10.1016/j.jmmm.2007.02.057) (2007). Copyright ©Else- 
vier B. V.l 



the peaks are centered at similar field values and they 
are quite narrow and almost symmetric around the min- 
imum. However, for the FC loops, apart for the obvious 
shift of the peak positions, the decreasing branch peak is 
symmetric and narrow, while the increasing branch peak 
is deeper and asymmetric, enclosing bigger area under 
the loop curve. 

Asymmetric loops are usually found in different bilay- 
ered syst ems^^^o^o^o^o7^m an d are also evi- 
dent in some core/shell nanoparticle systems. However, 
clear-cut experiments revealing the microscopic origin of 
this asymmetry have only been performed in the former 

case 127 1 149 1 209 1 210 Mogt theories f EB for thin filmS; al _ 

though considering the possibility of formation of domain 
walls during the magnetization reversal, are not able to 
account for origin of this asymmetry. Only in recent mi- 
cromagentic simulation o 182 ' 201 , an asymmetry has been 
observed. Also MC simulations of the DS model for a 
single^ or twined anisotropy axesi^ have shown that 
the observed asymmetries depend on the angle between 
the easy axis of the AFM and the applied magnetic field. 
More recently, hysteresis loops computed by MC simu- 
lations of a FM cylindrical dot in contact with an AFM 
based in a ferromagnetic domain wall model for the inter- 
facial coupling, exhibited also an asymmetric profile^. 

However, within the context of our model, in core/shell 
nanoparticles, the observed loop asymmetries arises 
solely by the competition between the interfacial ex- 
change coupling and the aligning effect of the magnetic 



field due to the intricate geometry at the interface. 



F. Reversal mechanisms 

These observations also indicate that the loop asym- 
metry reflects different reversal mechanisms in both 
branches of the hysteresis loops. This can be corrobo- 
rated by direct inspection of the spin configurations along 
the loops. In Fig. \7\ several snapshots of a midplane 
cross section parallel (left panel) and perpendicular (right 
panel) to the applied field direction, taken around the co- 
ercive fields Iiq for Ji nt = — 0.5 Jc are shown. As it is ev- 
idenced by the upper sequence of snapshots, the reversal 
proceeds by quasi uniform rotation along the descending 
branch, while nucleation of reversed domains at the inter- 
face and their subsequent propagation through the core 
center is basically the reversal process along the ascend- 
ing branch, as evidenced by the lower sequence in Fig. 
[3 Similar asymmetry between the loop branches has 
been also observed experimentally in bilayer a 127 i 149 i 209 
and, more recently, the relevance of nonuniform reversal 
modes to asymmetric magnetization reversal has been ev- 
idenced by measurements of hysteresis loops with vary- 
ing angle of the cooling field in Ni/NiO polycrystalline 
system 2 ^. A detailed inspection of the configurations, 
also reveal the presence of spins at core/shell interface 
aligned perpendicular to the field direction for inter- 
mediate field values (see for example the snapshots for 
h = —2.4,0.6 in Fig. [7J. This observation corroborates 
the interpretation of recent results of small-angle neu- 
tron scattering experiments on Fe oxidized nanoparticles, 
in which the anisotropy of the obtained spectra was at- 
tributed to the existence of a net magnetic component 
aligned perpendicularly to the field directio n 212 ' 213 . Note 
that similar perpendicular couplings have been observed 
in thin film system a 214 i 215 . 

The microscopic origin for the different reversal mech- 
anisms can be further clarified by looking at the behav- 
ior of the interface shell spins along the hysteresis loop 
(see Fig. 2}d,c). While in the descending branch there 
is a considerable amount of unpinned spins that are able 
to reverse following the core reversal, in the ascending 
branch remains constant (for Ji nt < 0), an indi- 

cation that spins at the shell interface remain pinned, 
hindering uniform rotation of the core but acting as a 
seed for the nucleation of reversed domains. 

The changes in the magnetic order at the core/shell 
interface and the presence of domain walls during reversal 
can be traced by monitoring the value of the average sum 
of the projection of the spin direction into the direction of 
the total magnetization vector along the hysteresis loops 
computed as 



m 1 



1 N 

(*) = jvEw-^w 



(6) 



This quantity should be close to 1 if the magnetization 
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FIG. 7: (Color online) Snapshots of the spin configurations of a midplane cross section of the particle parallel (upper sequence) 
and perpendicular (lower sequence) to the applied field direction taken at selected values of fields along the descending and 
ascending branches close to the coercive fields (/i^, /ij) for the case Ji n t = — 0.5Jc shown in Fig. la. 
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reversal proceeds by uniform rotation of the spins, since 
in this case the spins remain parallel to the global magne- 
tization direction. Deviations of m p (h) from 1 indicate 
the formation of non-uniform structures during the re- 
versal process. An example of the field variation of m p 
computed for all the core spins is shown in Fig. [BH, while 
in Fig. [8b we show m p (h) for the interfacial spins, where 
we have plotted separately the contribution of the core 
spins. 

During the decreasing field branch of the loop, m p re- 
mains quite close to 1 for the core spins, except for moder- 
ate decrease down for values of h close to . The sharp- 
ness and symmetry of the peak around confirms that 
the reversal proceeds by uniform rotation. In contrast, 
during the increasing field branch, an increasingly strong 
departure of m p from 1 starting from negative field values 
can be clearly observed, reaching its maximum value also 
near Kq. In this case, the observed peak asymmetry is 
indicative of the nucleation of the non-uniform domains 
observed in the snapshots of Fig. [7j These domains are 
formed at those points of the core interface with weaker 
values of the local exchange fields, as indicated by the 
more pronounced departure from 1 of m l p nt (h) (see Fig. 
[8b), than that corresponding to the total core magneti- 
zation (see Fig. [8^)- The variation of m p nt for interface 
shell spins during the decreasing branch indicates the ex- 
istence of a fraction of shell spins that reverse dragged 
by the spins at the core, while constancy of m p in the 
ascending branch is indicative of spins pinned during the 
core reversal. 

The origin of the loop asymmetry can be further clar- 
ified by monitoring the values the so-called overlap q(h) 
and link overlap qh{h) functions along the hysteresis 
loops, that are a generalization of similar quantities com- 
monly used in the spin-glass literatur e 216 ' 217 and that are 
defined as 




h (K) h (K) 



FIG. 8: (Color online) (Color online) Panel (a) shows the 
field dependence of the average spin projection of the core 
spins into the total magnetization direction Trip (squares). 
In panel (b), only the contribution of all the interface spins 
(circles) has been taken into account, while the contributions 
of core and shell spins at the interface are shown in squares 
and triangles, respectively. The field dependence of the link 
overlap §l and overlap q functions for the interfacial spins 
at the shell (circles) and at the core (squares) is shown in 
the panels at the right [(c) and (d), respectively], for Ji nt = 
-0.5 J c . 

of reversed nuclei at the particle core that sweep the par- 
ticle during reversal. 

The function q(h) measures differences of the spin con- 
figuration at field h with respect to the one attained after 
FC. Therefore, the decrease of q for the interface shell 
spins when reducing the magnetic field indicates the ex- 
istence of a fraction of shell spins that reverse dragged 
by core spins, while the constancy of q in the ascending 
branch reveals the existence of spins pinned during core 
reversal. 



i=l 

9* CO = E TT ' S 3 (h FC ) S t (h) ■ ^(h) , (7) 

(ij) 

where in qh(h) the summation is over nearest neighbors 
and Ni is a normalization factor that counts the number 
of bonds. 

An example of the field dependence of these overlaps, 
computed only for the interfacial spins, is shown in Figs. 
Hfc, d, where we have separated the contribution of the 
shell and core spins. A departure of q^ from unity is 
known to be proportional to the surface of reversed do- 
mains formed at field h and, therefore, gz, is sensible to 
the existence of non-uniform structures. The sharp de- 
crease of qh for core spins and the symmetry of the peak 
around the negative coercive field indicates uniform re- 
versal. However, the progressive reduction of (?l along 
the ascending branch and the asymmetry of the peak 
around the positive coercive field indicates the formation 



G. Vertical loop shifts 

Clearly correlated to the observation of EB and the 
loop asymmetry, the loops experience a shift along the 
vertical (M z ) axis which increases with Ji n t, as reflected 
in Fig. [6] by the difference of the M z values in the high 
field region of the two loop branches or at the remanence 
points. The field dependence of magnetization compo- 
nent transverse to the field direction, M tr (circles in the 
lower panels of Fig. [6]) , indicates that M tr attains values 
for the descending loop branch that are higher than in the 
ascending branch. Moreover, the M tr values around the 
peaks increase with increasing Ji„t . Snapshots of the spin 
configurations at the remanence points of the hysteresis 
loops are displayed in Fig. [9l They show the existence of 
a higher amount of core spins with transverse orientation 
near the interface at the lower branch (Fig. [9b, d) than at 
the upper branch. This is in agreement with the results 
of some experiments in oxidized particles where this ver- 
tical shift was also observeo^i^iZ 6 -^^ and with the 
observation of transverse magnetization components dur- 
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FIG. 9: (Color online) Snapshots of the remanent spin con- 
figurations of the upper (a,c panels) and lower (b,d panels) 
branches of the hystersis loops showing midplane cross sec- 
tions of the particle parallel (a,b) and perpendicular to the z 
axis (c, d) for the case Ji n t = — Jc shown in Fig. 1. Dark 
(light) blue cones represent core (core/interface) spins while 
green (yellow) ones are for spins at the shell (shell/interface). 



ing reversal revealed unambiguously by magneto-optical 
Kerr effect in bilayerai^. Our simulation results above, 
indicate that the microscopic origin of the vertical shift 
resides in the different reversal mechanisms on the two 
loop branches due to the existence of uncompensated 
pinned moments at the core/shell interface that facili- 
tate the nucleation of non-uniform magnetic structures 
during the ascending field branch of the loops. The re- 
cent observation that the vertical shift may be attributed 
to the existence defect moments^ will be checked within 
the scope of our model by removing some spins at the 
interface or at the core of the particle, this work is in 
progress. 




FIG. 10: (Color online) Thermal dependence of the coercive 
fields Hq , Kq, he and the exchange bias field h eb for a particle 
with Jint = — 0.5Jc. 



found a linear decrease of h e b vanishing below Tjv, in ex- 
cellent agreement with experimetal results^. 

The results of our finite temperature simulations for 
the particle with Ji n t = — 0.5Jc are displayed in Fig. 
lOl where the variation of , h^, he and /i b with the 
temperature at the end of the FC process are displayed. 
Let us notice first the different dependencies of on T. 
Starting from the lowest temperature, both quantities 
first decrease up to 2 K aprox. However, after further 



increase in T, whereas h c is stable up to Tg = 6 K, 
Hq increases, reaching a maximum at the same T. As 
a consequence, we find that /i b vanishes at 6 K, while 
he presents a maximum at the same temperature. This 
seems to agree with Trohidou's results^ for some of their 
particle sizes. This enhancement of he at the block- 
ing temperature Tb where /i e b vanishes has also been 
reported for bilayered systems, but, to our knowledge, 
not for particles. 



Other studies 



H. Temperature dependence of h e t 

The thermal dependence of h e b and he can also 
be studied by the MC method. Using a model for 
an oxidized nanoparticle similar to ours, Trohidou and 
coworkers^ first computed the thermal dependence 
of the coercive field founding that, compared to the 
non-oxidized particles, there was an increase of he in 
all the temperature range and also a reversal in the 
size dependence of the coercivity. They also found a 
steeper temperature dependence of he when the inter- 
face anisotropy is enhanced. More recently, they have 
also computed thermal dependences of /i c b for several 
particle size o 125 ' 185 ' 219 , finding a stronger temperature 
dependence for the bias field than for the coercive field. 
MC simulations of the DS model for bilayers^i have also 



The particle size and shell thickness dependence of EB 
has been studied by Trohidou and coworkers 1 - 8 ™ These 
authors argue that the observation of EB depends mainly 
on the structure of the interface and not on its size, also 
in agreement with our findings. They have found that a 
reduction of the core size for a given particle size enhances 
/icb and reduces he- The same group 74 ' 125 has recently 
performed simulations of a core/shell nanoparticle with 
random anisotropy directions in a FIM shell which repro- 
duce the experimentally found training effects in Fe oxi- 
dize particle o 70 ' 73 ' 126 and also the aging of the remanent 
magnetization, he and h e b- This last quantity was found 
to increse with the time during which the cooling field has 
been applied. The existence of training effects in layered 
systems has also been shown to be in agreement with ex- 
periments in MC simulations within the context of the 
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DS mode l 130 ' 181 . The role of imperfect interfaces in es- 
tablishing the EB has also been studied in this modeU^. 
Let us also mention the works by Usov at al . 221 ' 222 ' 223 , 
in which the authors compute the magnetic states and 
hysteresis loops of composite nanoparticles and bilayers 
using a quantum mechanical Hartree-Fock approxima- 
tion. Also using a quantum mechanical approach, Mata 
et alM^ have suggested that quantum zero temperature 
fluctuations of surface spins near an AFM surface induce 
dipole fields that may account for the observed exchange 
anisotropics. 

A microscopic model for interface roughness in bilay- 
ers was proposed by Almeida and RezendeA^, who com- 
puted the hysteresis loops for Ising spins in a mean-field 
approximation. Apart from the loop shift and enhanced 
coercivity, they showed that the sign of the exchange bias 
field changes as the initial temperature of the FC process 
is lowered and as the cooling field is varied, in agreement 
with experimental report o 136 ' 225 . In a model of bilayers 
based on a generalization of MB model that included bi- 
quadratic exchange and that accounted for the granular 
structure of the AFM, Hu et al . 226 ' 227 ' 228 computed the 
thermal dependence of he and h e h, in agreement with 
some experimental results. 

MC simulations by Lederman et al.^2 3 - for F9/F9F2 
demonstrated that EB is generated when the AF sub- 
lattices have an unequal exchange coupling with the FM 
and that perpendicular order between the FM and AFM 
is possible for large interface exchange coupling, in agree- 
ment with previous theorie o 138 ' 229 . In similar MC simu- 
lationsy, Billoni et aiM& have studied the influence of the 
value of exchange coupling constant at the interface on 
H e b at different temperatures. The effect of interfacial 
coupling on the magnetic ordering of models of coupled 
bilayers was studied using MC simulations by Tsai et 
aJ . 231 ' 232 and Alonso et a].— , by Finazzi 2 ^. in a micro- 
magnetic approach. 

Within the context of a random field Ising model, Ilia 
et a h 235 i 236 performed MC calculations of bilayers where 
the existence of EB was related to the fraction of en- 
hanced broken links and was shown to be due to mi- 
nor loop effects. The same approach has been used in a 
model that includes a partially covering of the FM/AFM 
interface by a non-magnetic spacer, showing its influence 
in perpendicular E B 237 ' 238 ' 239 . Also based on the same 
model, Meilikhov and Farzetdinova 2 ^ presented a mean- 
field approach that allows analytical solutions. MC sim- 
ulations of the related random anisotropy Ising model 
by Negulescu et alM^- showed also EB effects due to the 
roughness of the interface. 

First principle studies specifically addressing the origin 
of the EB effect are scarce. However, interesting calcu- 
lations of Co/FeMn bilayers by Nakamura et a\M2- using 
FLAPW method to incorporate noncollinear magnetic 
structures have demonstrated from first principles that 
an out-of-plane magnetic anisotropy is induced at the 



Co/FeMn interface, in accordance with experimental re- 
ports. 

V. CONCLUDING REMARKS 

We have reviewed the main phenomenology associated 
to EB in core/shell nanoparticles and presented details 
of our simulations of a model for these systems which 
explicitly takes into account the microscopic parameters 
characterizing the core/shell interface. The results of the 
simulations are able to account for some of the experi- 
mental observations. The obtained hysteresis loops after 
FC present shifts along the field axis which are directly 
related to the existence of a fraction of uncompensated 
spins at the shell interface that remain pinned during 
field cycling. The results of the simulations have revealed 
asymmetries in the hysteresis loops which, by detailed 
analysis of the microscopic magnetic configurations, have 
been linked to the occurrence of different reversal mech- 
anisms in the two loop branches. The existence of in- 
terfacial groups of spins aligned transverse to the field 
direction and the above mentioned difference in the re- 
versal mechanisms is also responsible for the vertical shift 
of the loops. Moreover, we have been able to establish 
a quantitative connection between the macroscopic mag- 
nitude of the EB fields and the microscopic value of the 
net magnetization of the interfacial shell spins. 

In order to account for the effects that other charac- 
teristic features of real nanoparticle systems may have 
on the experimentally observed phenomenology, further 
ingredients will have to be considered in simulations of 
microscopic models. Among them, let us mention the 
intrinsic surface spin disorder and surface anisotropy, 
the distribution in particle sizes and randomness of the 
anisotropy directions and the existence of dipolar inter- 
particle interactions in self-assembled or agglomerated 
particle assemblies. Finally, we hoped that the possi- 
bility to perform ab initio calculations of nanoscale clus- 
ters form first principles will lead to more realistic inputs 
for the microsocopic parameters needed for MC and mi- 
cromagnetics simulations, allowing a multiscale approach 
that will shed new light into the microscopic origin of the 
peculiar magnetic properties of nanoparticles. 
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